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Abstract 

Linearized complex Ginzburg-Landau equation models various physical phenomena 
and the stability controls of them are important. In this paper, we study the control of 
the LCGLE with a simultaneously space and time dependent coefficient by transforming 
it into a complex heat equation. It is shown that under certain conditions on the 
coefficient functions U2(x, t), the exponential stability of the system at any rate can be 
achieved by boundary control based on the state feedback. The kernels are explicitly 
calculated as series of approximation and shown to be twice differentiable by using the 
method of dominant. Both the exponential stabilities of the systems with Dirichlet and 
Neumann boundary conditions are strictly proven. 

Keywords: partial differential equations, Ginzburg-Landau equation, heat equation, bound- 
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1 Introduction 



In this paper, we use the boundary feedback control law, which is based on the back step- 
ping methodology, to stabilize the linearized complex Ginzburg-Landau equation (LCGLE) 

du(x.t) d 2 u(x,t) ..du(x,t) 

— ^ = ai-^P + a ^^— + a ^ *) ( L1 ) 

where x G (0,L), te (0,T), u : (0,L)x(0,T) -»• C, a x G C and &(ai) > 0, a 3 G ^([O^C) 
and a 2 G C 1 ([0, L] x [0, T];C) is a simultaneously space and time dependent coefficient 
function. The Dirichlet boundary conditions of the system are 

u(0,t) = p(t), (1.2) 
u{L,t) = (1.3) 

where p(t) is the control input. This LCGLE models various physical phenomena, such as 
the amplitude equation in pattern formation |2j and the reaction diffusion of two chemicals 
in one dimension |2JIH], in all of which the controls of the systems are important. 

Without losing any generality, we set a 3 (x) = since this term could always be eliminated 
by a gauge transform 

u -> u = ue- l,2 ^ a3{y)laidy (1.4) 
to obtain the equation we will concentrate on hereafter 

du(x,t) d 2 u(x,t) 

y =ai +a2{x,t)u{x,t). (1.5) 

Equation ()1.5j) is equivalent to a complex linear heat equation. The feedback control of this 
equation where the a\ and a 2 are real constants was first addressed by Boskovic, Krstic and 
Liu in 0] and the instability of the systems was also shown there. Later the plant coefficient 
was generalized to a space dependent case a 2 = a 2 (5;) by Liu |5.. More recently, Aamo, 
Smyshlyaev and Krstic PQ El El considered the generalization to a complex equation with a 
complex space or time dependent coefficient function, i.e., a 2 = a 2 (x) or a 2 = a 2 (t). In this 
paper, we extend their work by generalizing the plant coefficient to a simultaneously space 
and time dependent case, i.e., a 2 = a 2 (x, t). We also considered the Neumann boundary 
condition problems. Besides, from the proof of lemma (j2.1j) we can see that the explicit 
result for the control kernel is given and therefore it could be immediately used in numerical 
simulation. 

The paper is organized as follows. In sectional we prove that for some a 2 (x,t), the 
system (11.5(1 with the Dirichlet boundary conditions could be stabilized by the boundary 
feedback control. The unique existence of the kernels are shown in subsection 12.11 Using 
these kernels, in subsection 12.21 the original system is transformed into a new well-posed 
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system. In subsection 12.31 we show that the original system is exponentially stabilized. In 
section El we prove that the Neumann boundary problem can also be stabilized using the 
similar methodology. 



2 Dirichlet Boundary Value Problem 

First we break equation (jl.5|) into two coupled PDEs with real domains, real coefficient 
functions and real codomains by defining 



p(x,t) = U(u(x, t)) — - (u(x, t) + u(x, t)) , 
t(x,t) = Q(u(x, t)) — — (u(x, t) — u(x, t)) 

where 

x = — - — , t — — and u(x, t) = u(x, t) 
and ~ denotes the complex conjugate. Equation ()1.5|) then becomes 

Pt = clrPxx + b R {x, t)p - ajL xx - bj(x, t)i, 
k = a R L xx + b R (x, t)i + ajp xx + t)p 

for (x,t) G (0, 1) x (0, 1) = E, with boundary conditions 

p{0,t) = 0, 6(0, t) = 
p(l,t) = p R (t), t(l,t)= Pl (t) 



where 



a R = ai = j2^{ a i), 

b R = $t(a 2 (x,t)), hi = S(a 2 (5, t)). 



2.1 
2.2 

2.3 



2.4 
2.5 



2.6 
2.7 



2.8 
2.9 



Note that the transformation (J2.3j) is to normalize the x domain and the t domain. To 
stabilize the new equation system (j2.4|) - (j2.9J) . we set the boundary control input pit) into the 
form 

Pn(t) = [ [k(l,y,t)p(y,t) + k c (l,y,t)i(y,t)}dy (2.10 
Jo 

Pi(t) = [ [k{l,y,t)i{y,t)-k c {l,y,t)p{y,t)}dy, (2.11 
Jo 

where k and k r are kernels that we should find. 
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2.1 The Unique Existence of the Kernels 

First, we state a lemma about the kernels k and k c for future use. It will be clear in lemma 
(12. 2j) why we consider this PDE system. 

Lemma 2.1. Ifa R b R (x, t) + aibi(x, t) is analytic in t, the partial differential equation system 
of k and k c : 

kxx ^yy 

+ /3(x, y, i)k + (3 c (x, y, t)k c + p x k t + p 2 k c ,t (2.12) 

k c ,xx = k c>yy - Pc(x,y,t)k + P(x,y,t)k c + q 1 k c j + q 2 k t (2.13) 

for (x, y) G Q with boundary conditions 

k(x,x,t) = ~J (3(%%t)dry (2.14) 

k c (x,x,t) = /3 c ( 7 ,7,*)d 7 (2.15) 

k(x,0,t) = (2.16) 

k c (x,0,t) = (2.17) 

where 

Pi = <?i = a R /(a 2 R + a]), (2.18) 

P2 = -q 2 = -a I /{a 2 R + a 2 T ), (2.19) 

P(x,y,t) = M6 Ji ( 3 /,t)-/fl(x J t)) + ai(6j(3/,t)-/j(x,t))]/(a| l + a?), (2.20) 

(3 c {x,y,t) = [a R {b I {y,t)-f I (x,t))-a I (b R {y,t)-f R {x,t))}/(a 2 R + a 2 ), (2.21) 

m which f R (x,t), fi(y,t) defined in theorem \2.1)) are analytic in t and differentiate in y, 
has a unique solution satisfying 

\k{x,y,t)\ < Me M{x2 ~ y2) (2.22) 

\k c {x,y,t)\ < Me M ^~^ , (2.23) 

where M is a positive constant, for any time interval (0,to)(^o < 1) we are concerning. 
Proof. Using the substitutions 

x = £ + v, y = t~V, (2-24) 

G(£, T], t) = k(x, y, t), G c (£, ri, t) = k c (x, y, t), (2.25) 

5(f , rj, t) = P(x, y, t), 5 c (e, V, = y, *), ( 2 -26) 
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we can have 



k xx = 1/ 4((j££ + 2G^ V + G m ) k CtXX = 1/ 4(G Cj ^ + 2G C ^ V + G cm ) 

kyy — 1/4 (G^ — 2G^ V + G m ) k C}V y = l/4(G Cj ^ — 2G c ^ r) + G cm ) (2.27) 

kt = Gt k c j = G Cj t- 

Thus the system fl2.12j) - ()2.17|) is transformed into 

G^,V,t) = ^,r ] ,t)+p 1 ^G(^ V ,t)+^, V ,t)+p 2 ^G c ^ V ,^ (2.28) 
Gc£ V {Z,V,t) = (s(£,V,t) + qi^)G c (Z, V ,t)+ (-5 c (Z,rj,t) + q 2 ^)G(Z,ri,t) (2.29) 



with boundary conditions 



G(£,0,t) = -i / <)(t.(). t),)-. (2.30) 







G c (£,0,t) = \ j 5 c (r,0,t)6r } (2.31) 

G(£,£,t) = 0, (2.32) 

G e (£,t,t) = 0. (2.33) 

By direct integrating of (|2.28J) and (|2.29J) with respect to £ and rj and using (|2.30jl - (|2.31jl . we 
obtain the equivalent integral equations 

G{£, r),t) = ~J 8(t, 0, t)dr + y jT" (V, s, t) + Pl G(r, s, *)dsdr 

+ r r iw, *) + g ^ s > (2- 34 ) 



'77 </0 



G c (£, V,t) = \J S(t, 0, t)t*r + J jT (5(t, s, t) + qi ^j G c (t, s, t)dsdr 



+ J J* (-S c {r, s, t) + g 2 -| ) G(t, s, t)dsdr (2.35) 



To proceed, we set 



G (C, v, t) = -~ J S(t, 0, t)5r, G C)0 (£, 77, f) = ~ J S(r, 0, t)5r (2.36) 

G n+1 (£,rj,t) = J J ^6(T,s,t)+p r ^jG n (T,s,t)dsdr 

+ t T (s c (t, s, t) + p 2 ^\ G c , n (r, s, t)dsdr (2.37) 



G C)n+1 (£,r),t) = J j ^5(r,s,t) +Qi^j G c , n (r,s,t)dsdT 



'rj JO 

We want to show that both the series 



+ / r (-Sc{t, s, t) + g 2 -| ) G n (r, s, t)dsdr. (2.38) 




G = °n(^ V, t) and G C = J2 G c,n(^ V, t) (2.39) 



n=l n=l 

are uniquely convergent and thus by the method of successive approximation they are the 
solutions of system (l2~34]) - (l2~35T) . 

This could be done by the virtue of Colton's method of dominant (HI E]. This method 
works as follows. If we are given two series 

oo oo 

Si = J> ln r, S 2 = J2 a 2nt n , t G (0, 1), (2.40) 

n=l n=l 

where a 2n > 0, then we say S 2 dominates Si if |oi n | < am-, n — 1, 2, 3, • • • , and write Si <C S 2 . 
It can be easily checked that 

ifSi<S 2 , then|Si|<S 2 , (2.41) 

and « S 1 « S 2 (l - ty l ; (2.42) 

if Sx < S 2 , S 2 < S 3 , then S x < S 3 ; (2.43) 
if S x < S 2 , S 3 < S 4 , then Sx + S 2 < S 3 + S 4 . (2.44) 

Moreover, the property of dominancy can also be kept if the integrals of the two series are 
not with respect to t but other variables, that is 

if Si < S 2 , then / Sidx < / S 2 dx. (2.45) 

On the other hand, if a function / is analytic in t 6 (0, 1), then there exist a positive constant 
C such that / < C(l - t)" 1 . 

Using this method in our problem, it can be shown that equations (|2.12J) and (|2.13|) have 
unique twice continuously differentiable solutions if S and S c are analytic in t. In fact, since 
o-rOr^x, t) +aibj(x, t) is analytic in t and fi(x, t) and Jr(x, t) can also be chosen to be analytic 
in t, from equation ()2.20|) . ()2.21|) and ()2.26|) . we know that 8 and 5 C are analytic in t. Thus 
we can let C, C c be two positive constants larger than 1 such that 

5(e,r ? ,t)«C(l-t)- 1 , 5 c (£,^)«C c (l-t)-\ (2.46) 
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or further 

5(£, v, t) < N{1 - ty\ 5 c (t rj, t) < N(l - t)~\ N = max{C, C c }. (2.47) 

From equation of ()2.36|) and noticing ()2.45|) . we have 

G < Nil-t)- 1 (2.48) 
G cfl < Nil-t)- 1 . (2.49) 

By induction, suppose 

G n « ^T Nn+ \ l _ t y{n + i) (25Q) 
n! 

G c , n « (2-51) 
n! 

we will have from ()2.37|) 

G n+ i= J J ^6 + pt^j G n dsdr + J J ^5 C + P2^j G^dsdr 

« „ ^ nn2 (1 - t)- (n+2 \2N + npi + np 2 ) 

« 4 " +i f:n;r 2 ' 1 - ^ 



where all the properties ()2.41|) - ()2.45|) have been used. Similarly, we also have 

[n + 1)! 

and hence by ()2.41|) 



|GWrl < 4 H-r# (1 " t)- ( " +2) , |<WI < 4 ™ * (1 - t)"^. (2.54) 
(n + 1)! (n + 1)! 

It is clear that the two series convergent absolutely and uniformly and are solution of 
fl2~Ml) - (l2~35J) . G and G c are C 2 since S(£,r],t) and S c (£,r),t) are C 1 . Since 

/ N \ n 

A n {^YN n {\ - t)- n = Ux 2 - y 2 )— t ) > ( 2 - 55 ) 

(I2~2"2l and (l2~2~31) follow directly by assigning N/(l - t Q ) = M. □ 
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Remark 2.1. (1). The proof of Lemma X2.1)) provides a numerical computation scheme of 
successive approximation to compute the kernels. This makes the feedback laws very useful 
in real problems. (2). We require t G (0,t ), t < 1 in \2. 22]) and \2.23j) since from \2.5Sty 
one can see that k and k c are not bounded when t — > 1. This requirement is tolerable since in 
practice, we only require the system to be stable in a time interval t G (0, To) where Tq = t^T 
from rescaling equation \2. , 6 j\) and the T can be choose as large as we want. 

2.2 Conversion and Inverse Convertibility of the Original System 

We want to show that the original system (|2.4j) - (|2.11j) can be converted to a new system by 
integral transform with k and k c as kernels and this conversion is invertible. 

Lemma 2.2. Let k(x,y,t) and k c (x,y,t) be the solution of A2.12)) - f2.1ty) and define a pair 
of linear bounded operator K and K c by 



p(x,t) = (Kp)(x,t) = p(x,t) - [k(x,y,t)p(y,t) + k c (x,y,t)i(y,t)]dy (2.56) 

Jo 

/•x 

l(x,t) = (K c i)(x,t) = i(x,t) - [-k c (x,y,t)p(y,t) + k(x,y,t)i(y,t)]dy, (2.57) 



o 



Then, 

1. K and K c convert system \2.4]) - ^Tl]) to system 

Pt = cirPxx + fii(x, t)p - ail xx - fi(x, t)l, (2.58) 
k = cbipxx + fi(x, t)p + a R l xx + f R (x, t)Z, (2.59) 

for x G (0,1), with boundary conditions 

p(0,t) = 0, 1(0, t) = 0, p(l,t) = 0, 1(1, t) = 0. (2.60) 

2. Both K and K c have linear bounded inverses. 

Proof. Proof of part 1.: To prove 1, we compute as follows. Differentiating ()2.56|) and using 
Q-flZinj), we have 

p t (x, t) =a R p xx + b R (x, t)p - ajL xx - bj(x, t)i 
f ' 

- / [k(x,y,t)(a R p xx + b R (y,t)p-a I L xx -b I (y,t)L) + k t p 
Jo 

k c (x,y,t)(a!p xx + bj(y, t)p + o r l xx + b R (y,t)i) + k c ^]dy. 
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Using (|2.56|) - (j2.57|) . integrating the integral by parts twice and using (|2.6|) . (j2.16j) - (|2.17|) and 
(ETEUl) . we will get 



,2 r x 



p t (x, t) =a R (p xx (x, t) + ^ J [k(x, y, t)p(y, t) + k c (x, y, t)t(y, t)]dy 
+b R (x, t) (p(x, t)+ J [k(x, y, t)p(y, t) + k c (x, y, t)i(y, t)]dy^j 

( d 2 r 

\L xx (x,t) + — J [-k c (x,y,t)p(y,t)+k(x,y,t)i(y,t))dy 
t) ( T(x, t)+ [~k c (x, y, t)p(y, t) + k(x, y, t)i{y, t)]dy 



o 



- k(x, x, t)a R p x (x, t) + k(x, x, t)aji x (x, t) - k c (x, x, t)ajp x (x, t) - k c (x, x, t)a R i x (x, t) 
+ k y (x, x, t)a R p(x, t) — k y (x, x, t)aji{x, t) + k c>y (x, x, t)ajp(x, t) + k c>y (x, x, t)a R i(x, t) 

[k yy (x, y, t)(a R p(y, t) - an(y, t)) + k(x, y, t)(b R (y, t)p(y, t) - bj(y, t)t(y, t)) 

+ k c>yy (x, y, t)(a/p(y, t) + a R i(y, t)) + k c (x, y, t)(6/(y, t)p(y, t) + b R (y, t)t(y, t)) 
+ h(x, y, t)p(y, t) + k Cjt (x, y, t)i(y, t)]dy. 



o 



(2.61) 



Using the identity 

2 



dx 2 



Kx, y, t)x(y, t)dy = / X xx (x, y, t)x(y, t)dy + X x (x, x, t)x(x, t) 

i! JO 

+ X(x, t) ^^^ + A(x, x, t) Xx (x, t) (2.62) 



and the equation (|2.58j) - (|2.59|) again, we obtain 

p t (x, t) =a R p xx (x, t) - afi xx {x, t) + b R (x, t)p(x, t) - 6/(x, t)l(x, t) 

+ 2 (a / fc( ^ M) + ai d J^lA ^ ~ p{Xj t) + £ R{X} y> t)p(Vi t)dy 

+ 2 / aR dk c (x x,t) _ a dk( x x ,t) \ + r ^ ^ t)L ^ (2 _ 63) 



\ dx dx j j 

where 

R(x, y, t) =a R (k xx (x, y, t) - k yy (x, y, t)) + ai{k c ^ xx {x, y, t) - k c , yy (x, y, £)) - k t (x, y, t) 

„/ dk(x,x,t) dk c (x,x,t) , , . , . A , , . 
+ 2U R dx + a i dx J +b R {x,t)-b R {y,t)j k(x,y,t) 

( dk c (x,x,t) dk(x,x,t) , . . , . A , / 
+ 2 ( -a R cK dx 1 + aj y Qx ! + bj(x, t) - b T {y, t) j k c (x, y, t), (2.64) 
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I(x, y, t) =a R (k CtXX (x, y, t) - k c , yy (x, y, £)) - ai{k xx {x, y, t) - k yy (x, y, t)) - k C)t (x, y, t) 
+ 2 ^ CLR dk ( x ^ t ) + a dk e (x,x,t) + b ^ t) _ ^ ^ ^ ^ t) 

/ dk c (x,x,t) dk(x,x,t) , . . , . A , , . 
+ 2 (a R cy dx ' - aj > - b T {x, t) + b T {y, t) k(x, y, t). 

(2.65) 

Substituting equation (J2~T2|) - (12~T3|) and $ZTM - $TEB into (ESI) -(EES), we can get 

R(x, y, t) = [(a R pi + aiq 2 ) - l]k t (x, y, t) + (a R p 2 + aiqi)k Ctt (x, y, £)) 

+ (a R /3(x, y, t) - a//3 c (x, y, t) - (a R (3(x, x, t) - ai/3 c (x, x, t)) + b R (x, t) - b R (y, t))k(x, y, t) 

+ {a R p c (x, y, t) + a I f3 c {x, y, t) - {a R /3 c (x, x, t) + ai/3(x, x, t)) + 6j(x, t) - 6/(y, t))k c (x, y, t), 

(2.66) 

I(x, y, t) = [(a R qi - am) - l)k C;t (x, y, t) + (a R q 2 + a^k^x, y, £)) 

+ {-a R /3 c (x, y, t) - arfix, y, t) + {a R (3 c (x, x, t) + a T p(x, x, t)) - 6/(x, t) + h{y, t))k(x, y, t) 

+ {a R (3(x, y, t) - aiP c (x, y, t) - (a R (3(x, x, t) - ai/3 e (x, x, t)) + b R (x, t) - b R (y, t))k c (x, y, t). 

(2.67) 

From (j2~27H) - (l2~2T[) . we have 

a R p(x,y,t) - aj(3 c (x,y,t) = b R (y,t) - f R (x,t) (2.68) 
a R (3 c (x,y,t) +aiP{x,y,t) = bi(y,t) - fi{x,t) (2.69) 

and with the help of (j2.18|) - (|2.19|) . it is easy to check that 

R(x,y,t) = l(x,y,t) = 0. (2.70) 

Thus equation (J2.63j) becomes the following after using (|2.14|) . (|2.15|) again, 

p(x,t) =a R p xx (x,t) - ail xx (x,t) 

+ (—a R /3(x, x, t) + aij3 c (x, x, t) + b R (x, t))p(x, t) 

+ ( a R/3c(%, x, t) + aj/3(x, x, t) — bj(x, t))Z(x, t). (2-71) 

Substituting (J2.68[) and ()2.69|) into above, we can obtain ([2.580 . Equation (j2.59[) follows 
similarly by differentiating ([2.57]) and applying the same procedure. The boundary conditions 
(j2.60jl follow from setting x = and x — 1 in (j2.56|) and (j2.57|) and using (j2.6|) - ([2.7|) and 
([2~TT1jl - ([2~TT]l . Thus we proved part 1. 
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Proof of part 2.: To prove 2, the invertibility of the transform, consider the following 
transform 

px 

p(x,t) = (K- l p)(x,t)=p{x,t)- [l{x,y,t)p(y,t) + l c {x,y,t)i{y,t)]dy (2.72) 

Jo 

/•x 

t(x,t) = (K~ 1 t)(x,t)=Z(x,t)- I [-l c (x,y,t)p(y,t) +l(x,y,t)i(y,t)]dy, . (2.73) 

Jo 

They are actually the inverse transform if there exist unique solutions for l(x,y,t) and 
l c (x,y,t) that satisfying 

Ixx = k yy ~ P(y,x,t)k - (3 c (y,x,t)k c - Pik t - p 2 k c ,t (2.74) 
l c ,xx = k Ct yy + (3 c (y,x,t)k - /3(y,x,t)k c - qik Cjt - q 2 k t (2.75) 

with boundary conditions 

l(x, x, t) 

Z c (x, x, t) 
l(x,0,t) 

and 

\l(x,y,t)\ < Me M{x2 - y2 \ (2.79) 
\l c (x,y,t)\ < Me M ^~^ , (2.80) 

where M is a positive constant, for any time interval (0, to)(^o < 1) we are concerning. 

The proof that the transform (|2.72j) and (|2.73|) will lead to system (|2.74)) - (|2.78)) is similar 
to the proof of part 1 of this lemma. It is also clear that the transforms are linear. The 
unique existence of the solution to system ([2.74p - ([2.8U|) is similar to the proof of lemma (j2.1j) 
and thus the transforms are also bounded. Thus this proved that there exist the bounded 
linear inverse transforms as shown in (J2.72)) and (J2.73j) . □ 

2.3 Stability of the Controlled System 

Our purpose of all doing this is to control the original system by the boundary inputs. Its 
stability is established in the following theorem. 

Theorem 2.1. There exist feedback kernels k(l,-,t), k c (l,-,t) E C 2,oo (0, 1) x (0,i ),*o < 1, 
such that for arbitrary initial data p°(x), t°(x) G ff 1 (0, 1), system l[2.4{ )- l[2. 7| ) with \2.11\) - 
\2.1(J\) . where bji(x,t) and bi(x,t) are analytic in t, has a unique solution p, i € C 2,oo (0, 1) x 
(0, to) that is exponentially stable in the L 2 (0, 1) and Hi(0, 1) norms. 



P c (j,j,t)dj 



2 jo 
Ux,0,t) = 



(2.76) 

(2.77) 
(2.78) 
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Proof. First we notice that problem (|2.4j) - (j2.7|) with (|2.11|) - (|2.1(J|) is well posed since, by 
lemma (j2.2j) they can be transformed into the problem ()2.58|) - (j2.60|) via the isomorphism 
defined by (l23oT) and (12371) and problem (I23HJ) - (12301) is well posed (see, e.g. [TUj). 

Since through part 2 of lemma ()2.2|) . we showed that there exist linear bounded invertible 
transforms, it is sufficient to show that (p, I) are exponentially stable. The proof of this could 
be done by the same definitions and evaluations of the energy and potential as in corollary 
7 of Ref. [I]. Especially the fj(x,t), f R (x,t) used in lemma ()2.1|) are also defined there. □ 

3 Neumann Boundary Value Problem 

To stabilize the Neumann boundary value problem 

Pt = a R p xx + b R (x,t)p- a I L xx -b I (x,t)L, (3.1 
k = a R i xx + b R (x,t)i + a I p xx + b I (x,t)p (3.2 

for x G (0, l),t G (0, 1), with boundary conditions 

p x {0,t) = 0, t x {0,t) = (3.3 
p(l,t) = p R (t), L(l,t)= Pl (t) , (3.4 

the transform as (|2.5fij) - (|2.57|) with kernel k(x,y,t) and k c (x,y,t) in the following system 

k xx = k y y + P(x,y,t)k + p c (x,y,t)k c + pik t +p 2 k Ctt (3.5 
k c ,xx = k c , yy - (3 c (x,y,t)k + P(x,y,t)k c + q 1 k Ctt + q 2 k t (3.6 



for (x, y) G fl with boundary conditions 

~2 7 



k(x,x,t) = ~ I /%,7,t)d7 (3.7 



1 



X 



k c (x,x,t) = - J /? c (7,7,t)c?7 (3.< l 



o 



k x (x,0,t) = k CiX {x,0,t) = (3.9 
k{0,0,t) = fc c (0,0,t) = (3.10 

where P, Pc,Pi,P2,Qi,Q2, fi and f R are the same as in lemma (j2.1j) . will lead to the same 
system as ()2.58|) - ()2.60|) . If we can show that the above kernels uniquely exist, all the other 
lemmas and theorem in section |21 can be applied in this case. The main result are established 
as follows, in which the unique existence of the new kernels are explicitly shown. 
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Lemma 3.1. If a^b^x, t) + ajbi(x, t) is analytic in t, the partial differential equation system 
\3.ty) - K3~Ttyl has a unique solution satisfying 



\k{x,y,t)\ < Me M(x2 - y2) 
\k c (x,y,t)\ < Me M ^- y2 \ 

where M is a positive constant for t G (0, to) we are concerning. 



(3.11) 
(3.12) 



Proof. The proof of this lemma is similar to that of lemma ()2.1|) except we should find 
a G(£,£,t) as ()2.32j) in lemma (|2.1j) since this must be used in order to obtain G(^,t],t) as 
(j2.34j) . Using the same substitutions as in ()2.24|) - (j2.26j) . system ()3.7|) - (j3.8j) can be transformed 
into 



G 



d 



G 



cfy 



dt 




d 



<K£, V,t)+Piiz) G(Z, V,t)+[ 5 C (£ } 77, t) + pa- G c (£, V, *) (3-13) 



dt 



with boundary conditions 



dt 



G{t,0,t) 
G c (Z,0,t) 

GcM,U) 
G(0,0,t) 







6(Z, rj,t)+ qi -\ G c {t, T), t) + -6 C (£, V: t) + q 2 -\ % t) (3.14) 



dt 



8(r,0,t)dr, 



6 c (r,0,t)dr, 

G c ,T){ii i-, t) 

G c (0,0,t) = 0. 



Differentiating equation ()3.15|) with respect to £ gives 

G c (£,0,t) = 0,t). 
Integrating equation (|3.13|) with respect to rj from to £ and using (|3.2(Jj) gives 

=G^, 0,t) + J o s, t) + Pl G{£, s, t) + (s c (£, s, t) + p 2 G c (£, s, t) 



(3.15) 

(3.16) 

(3.17) 
(3.18) 
(3.19) 



(3.20) 



ds 



i r 

-(5(e,o,t) + 



5 



at 



, 5, t) + pi^- G(£, s, *) + 5 C (£, s, t) + p 27 - G c (£, a, t) 



d 



■dt 



ds. 

(3.21) 
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Thus using (|3.17|) . we have 
dG(Z,Z,t) 

d 



-5(£,0,t)+2 



dt 



^t)+ Pl —) s,t)+[ 5 c (f , s, t) + p 2 — G c {£, s, t) 







■dt 



ds. 
(3.22) 

Integrating the above equation from to £ and using (j3.19|) gives the function we wanted 



6(r,0,t)dT 



+ 2 




o J o 







dt 



S(t, s,t)+p 1 —\ G(t, s, t) + 5 c (t, s, t) + p 2 — G c (t, s, t) 







'dt 



dsdr. (3.23) 



Using ()3.23|) and integrating twice ()3.13|) first with respect to rj from to rj and second with 
respect to £ from rj to £, we can have the following integral equation 



S(r,0,t)dr 



v r T 



+ 2 



+ 




o jo 







dt 



5(T, S,t)+Pi— G(T, S, t) + 5 c (t, S, t) + G c (t, S, t) 



d 



■dt 




rj JO 



d 



dt 



5(t, s,t)+p 1 — G(t, s,t) + [ 5 c (t, s, t) + ^2777 G c (t, s, t) 



d 



'dt 



dsdr 
dsdr. (3.24) 



Similarly, for G c (£,rj,t) we will have 



6(r,0,t)dT 



+ 2 



+ 




o Jo 



d 



dt 



5(r,s,t) + qiT^ G c (r,s,t) + -8 c (r,s,t) + ^277; G(T,s,t) 



d 



■dt 



dsdr 




n Jo 



d 



dt 



6(t, s,t) + q 1 —) G c (t, s,t)+[ S c (t, s, t) + q 21 - G(r, s, t) 



d 



'dt 



dsdr. (3.25) 



Similarly as in lemma (j2.1|) . the above integral equations can also be shown to have uniquely 
and uniformly convergent solution satisfying (j3.11)) and (|3.12jl . □ 

Lemma 3.2. Let k(x,y,t) and k c (x,y,t) be the solution of ([.*?. .5)) - ( TP)) and define a pair of 
linear bounded operator K and K c by 

r 

p(x,t) = p(x,t)- [k(x,y,t)p(y,t) + k c {x,y,t)i(y,t)]dy (3.26) 
Jo 

Z(x,t) = i(x,t) - [-k c (x,y,t)p(y,t)+k(x,y,t)t(y,t)]dy, (3.27) 
Jo 

14 



Then, 

1. k and k c convert system \3. l\t - \3.$ to system 

Pt = urPxx + Sr{x, t)p - afi xx - //(x, t)l, (3.28) 

It = ai~p xx + f T (x, t)p + a R l xx + f R (x, t)i, (3.29) 

for x G (0, 1), with boundary conditions 

p(0, t) = 0, Z(0, t) = 0, p(l, t) = 0, fc(l, t) = 0. (3.30) 

2. Both K and K c have linear bounded inverses. 

Proof. The proof is similar to the proof of lemma fl2.2(l . □ 

Theorem 3.1. There exist feedback kernel k(l,-,t), k c (l,-,t) G C 2 '°°(0, 1) x (0,t ),*o < !, 
siic/i t/iat /or arbitrary initial data p°(x), l°(x) G -H' 1 (0, 1), system \3. l\i - \3.4\l with \2.11\) - 

has a unique solution p, i G C 2 '°°(0, 1) x (0, to) ^ s exponentially stable in the 
L 2 (0, 1) and H 1 ^, 1) norms. 

Proof. The proof is similar to the proof of theorem (J2.1j) . □ 



4 Remarks 

A much more challenging but also more important task is to stabilize the nonlinear heat 
equation with a u(x, t) dependent coefficient 

du(x,t) d 2 u(x,t) ~ _ ~ _ ~ 

— ^ — =ai — ^2 \-a2{x,t,u{x,t))u{x,t), (4.1) 

which models virous diffusion processes with a self-dependent source for the diffused matter by 
using the boundary feedback control method. It could be interesting to study the applicability 
of the method of dominant in this case. 
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